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INTRODUCTION 


The  design  of  improved  squibs  and  detonators  depends  on  better  understanding  the 
problem  of  the  initiation  of  exothermic  materials  by  imbedded  hot  wires.  Existing 
thermal  explosion  theory  1/  has  been  shown  to  predict  experimental  hot-wire  ignition 
delay  time,  but  does  not  show  the  dependence  of  ignition  delay  on  the  form  of  heating 
function  and  the  choice  of  explosive  and  wire  material.  Therefore,  the  reactivity  of  the 
imbedding  medium  has  been  characterized  by  an  "ignition  temperature"  rather  than  by 
the  more  fundamental  kinetic  parameters  of  thermal  explosion  theory.  Works  2/  and  3/ 
based  on  the  concept  of  the  fixed  ignition  temperature  shows  that  ignition  temperature 
depends  on  the  experimental  conditions  as  well  as  the  chosen  wire-explosive  system.  The 
ignition  temperature  of  a  composite  propellant  heated  by  current  pulses  at  preset 
duration  has  been  investigated  4/  and  5/,  and  it  was  found  that  this  explosion  property 
depends  on  ignition  delay  and  wire  diameter. 

A  different  thermal  explosion  criterion  6/,  which  incorporates  the  parametric 
dependence  of  the  delay  time  on  the  form  of  heating  function  and  the  choice  of  explosive 
and  wire  material, has  been  developed.  In  developing  the  theoretical  work  for  this 
explosion  criterion,  the  system  modeled  was  considered  an  exothermic  material 
surrounding  an  infinitely  long  heated  wire  of  infinite  thermal  diffusivity  and  of  radius  a . 
The  thermal  explosion  criterion  assumes  that  the  nominal  explosion  site  is  the 
wire/explasive  interface  r  =  a. The  interface  temperature  was  an  approximation  to  the 
solution  of  a  long  cylinder  in  perfect  contact  with  perfect  conductor  (the  wire).  The 
explosion  temperature  Tewas  that  value  of  T  which  minimizes  the  function  F(T). 


F(T)  =  trm(T)  +  tad(T)  (1) 


1/  Mauger,  F.,  The  Theory  of  Thermal  Explosions;  The  Initiation  of  an  Explosion  bv  Hot 
Wire,  ARDE  Report  No.  9168,  1960. 

2/  Altman,  D.  and  Grant,  Jr.,  A.  F.,  Fourth  Symposium  (International)  on  Combustion, 
(Baltimore  Williams  and  Wilkins,  1953)  p.  158. 

3/  Kabin,  I.,  Rosenthal,  L.  A.  and  Salem,  A.  D.,  The  Response  of  Electro-Explosive 
Devices  to  Transient  Electrical  Pulses,  N  OL  TR  61-20,  1961. 

4/  Baer,  A.  D.,  "Pyrodynamies",  3,  15  (1965). 

5/  Cheng,  J.,  Bouch,  L.  S.,  Keller,  J.  A.,  Baer,  A.  D.  and  Ryan,  N.  W.,  Ignition  and 
Combustion  of  Solid  Propellants,  AFOSR  66-1672,  1965.  ~ 

6/  Friedman,  M.  H.,  "A  General  Thermal  Explosion  Criterion  Application  to  Initiation  by 
Imbedded  Wires."  Combustion  <5c  Flame,  December  1969. 


where  tm(T)  is  the  time  at  which  the  interface  temperature  is  T  ,  and  tid  is  the 
adiabatic  explosion  time  explained  by  Friedman  6/.  Contrary  to  practical  consideration, 
the  analysis  cited  above  did  not  consider  the  contact  resistance  between  the  wire  and  the 
explosive  in  calculating  the  interface  temperature.  It  is  expected  that  due  to  contact 
resistance,  the  interface  temperature  will  lag  behind  the  wire  temperature,  resulting  in 
longer  ignition  delays  than  previously  predicted. 

Blackwell  7/  calculated  the  wire  temperature  caused  by  heat  supplied  at  constant 
rate.  He  also  developed  a  solution  to  the  surrounding  medium  temperature  in  the  Laplace 
Transform  space.  These  solutions  are  in  terms  of  an  infinite  integral  containing 
complicated  expressions  of  Bessel  functions  of  the  first  and  the  second  kind.  This 
method  was  used  to  develop  a  more  general  solution  for  the  temperature-time  history  of 
the  interface.  This  solution  is  valid  for  variable  current  heating  as  well  as  for  constant 
current  heating.  It  incorporates  the  parametric  dependence  of  the  delay  time  on  the 
contact  resistance  between  the  wire  and  the  explosive.  This  solution  is  used  with  the 
explosion  criterion  to  more  accurately  predict  ignition  delay  times. 

THEORETICAL  ANALYSIS 

Assume  an  infinite  region  of  one  material  (the  explosive)  bounded  internally  by  a 
circular  cylinder  of  another  perfectly  conducting  material  (the  wire).  With  the  initial  and 
boundary  conditions, 

a)  Zero  initial  temperature  throughout 

b)  Variable  heat  input  per  unit  time  to  the  inner  cylinder 

c)  Contact-resistance  at  the  boundary  between  cylinder  and  external  medium,  then 


d2e2  i  m2  i  ae2 
~drr+  *  ~  hjf"*" 

Initial  Conditions 


b<r<»,  t>0 


(2) 


0|  ■  02  a  0  at  t  =  0 


(3) 


7/  Blackwell,  J.  H.,  "A  Transient-Flow  Method  for  Determination  of  Thermal  Constants 
”  of  Insulating  Materials  in  Bulk,”  Journal  of  Applied  Physics,  Vol.  25,  No.  2,  1954. 


Boundary  Conditions 
6  2  is  bounded  as  r  -*•  «> 


ae,  „  ae. 

-K!  aT'Q’o-flT 


r  =  b,  t  >o 


r  =  b,  t  >o 


where  dt(t)  and  02(t,r)  are  the  temperature  of  the  wire  and  the  explosive  at  any  time, 
and  radius  r  is  the  radial  coordinate,  h|  is  the  diffusivity  of  the  explosive  material,  K2 
is  the  thermal  conductivity  of  the  explosive  material,  b  is  the  radius  of  the  wire,  t  is  the 
time,  H  is  the  outer  conductivity  at  the  surface  r  *  b. 


a  *  Mj  Cj  /  irb 


Q'(t)  =  Q(t)/2trb 


where  Mj  and  Cj  are  the  mass/unit  length  and  the  specific  heat  of  the  wire,  andQ(t)is 
the  heat  supplied/unit  wire  length/unit  time. 


Taking  Laplace  Transform  of  the  differential  equation  and  the  boundary  conditions 
results  in: 


d202  j  d©2 
dr2  r  dr 


q202 


(5) 


~K2  H<0i-02> 


Q'(P)  -  y  0i 


r  =  b 


(7) 


where  e^P) ,  02(P)  are  the  Laplace  Transforms  of  5,(t)  and  d2 (t , r), respectively,  p 
is  the  transform  variable,  and  q  *  >/F  /  hj . 
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The  solution  of  equation  5,  subject  to  the  boundary  conditions  of  equations  6  and  7,  is 


e  m  =  _ 2  bQ'(P)A(b) 

lV  [PabA(b)  +  2K2K1(qb>] 


_ 2  bQ'(P)  K0(qr) 

©2(P,  r>  =  qb  ( Pab  A(b)  +  2  K2  Kj  (qb)  ] 


Kc(qr)  /K2\ 

A(r)  =  ~  +  /  Kl  CQr)  ’  K° (<lr) ’  Ki 


where  qb  \  bH/  1  0  1  are  modified  Bessel  functions  of  the 

second  kind,  and  zero  and  first  orders  respectively,  and  A(b)  is  A(r)  at  r  =  b. 
Substituting  r  =  b  into  equation  9  gives  the  Laplace  Transform  of  the  interface 
temperature,  02(P,b).  Q^Pjj  and  02(P)  can  be  written  as 

01(P)  =  F1(P)*Q'(P)  (10) 


02(P)  *  F2(P)  •  Q'(P) 


where 


F  (P)  »  _ ***» _ 

(PabA(b)  +  2K2K1(qb)] 


2  b  K0  (qb) 

F2(P)  "  qb  [PobA(b)  +  2  K2 Kj  (qb)] 


Applying  the  inversion  theorem  of  the  Laplace  Transformation,  we  get 


i  ry“" 

fi(t)  *  57i  /  F.,p)'*P,dP 

Jr-i~ 


6 


F2(P)e+p«dp 


(15) 


f2(t) 


The  integrands  of  equations  14  and  15,  each  has  a  branch  point  at  the  origin,  and  neither 
has  poles  within  or  on  the  closed  contour  AKJGFEDCBA  .  Denoting  this  contour  by  C  , 
we  may  write 


si  fF-<P)e',dp 


(16) 


(17) 


By  Cauchy’s  Theorem,  the  integrals  on  the  left  side  of  equations  16  and  17  are  zeros.  As 
the  radius  R  of  the  outside  circular  arc  -*■  o®  ,  and  as  the  radius  e  of  the  inside  circular 
arc  -*>  0  *  the  second,  the  fourth,  and  the  sixth  integrals  of  the  right  side  of  each, 
equations  16  and  17  vanish.  Thus 


(24) 


Pj  =  xJ0(bx)  +  K2  Jx  (bx)  |  ^  x2  -  — I 
Qt  =  x  Y0(bx)  +  K2  Yj  (bx)  J-^  x2  -  -~J 


(25) 


where  J0(bx),  Jt(bx)  are  Bessel  functions  of  first  kind  of  zero  and  first  order 
respectively,  and  Y0(bx),  Yt(bx)  are  Bessel  functions  of  the  second  kind  of  zero  and  first 
order,  respectively.  Equations  10  and  11  show  that  the  Laplace  Transform  of  the  wire 
temperature  and  the  interface  temperature  are  the  product  of  two  transforms 
each  Fj  (P)  •  Q'(P)  ,  and  F2(P)  •  Q'(P)  .respectively.  A  convolution  solution  is  possible 

in  principle. 


0j(t)  = 


02(t)  = 


f^iOQ'a-  u)  du 


f2(u)  Q'(t-  u)  du 


substituting  equations  22  and  23  into  equations  26  and  27,  we  get 


(26) 


(27) 


0j(t)  = 


16  K2 
a2b  ir2h^ 


e  X  hzU  Q'(t-  u)  dx  du 
x  [  Pf  +  Qj  ] 


(28) 


02(t,b) 


16  K2 
a2  b  jt2 


«h2*2\  -*2h2u 

2H _L__ 


Q'(t  -  u)  dx 


x[P2  +  Q2] 


du 


(29) 
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The  solution  for  constant  current  heating  may  be  obtained  by  replacing  Q'(t)  in 
equations  28  and  29  by  the  constant  value  Q' ,  and  perform  the  integration  with  respect  to 
the  dummy  variable  u  . 


,ro.2!52L  /*•<'-«-**> 

1  a2  b  7r2  h,  -l  x3  [  P2  +  Qj2  ] 


02(t,b)  = 


/  al 


-x*hi t 


a2  b  7r2  h* 


x3  [Pi  +Qil 


Equation  30  is  the  solution  derived  by  Blackwell  7/  for  the  probe  temperature.  For  zero 
contact-resistance  between  the  wire  and  the  explosive,  the  solution  may  be  obtained  by 
letting  H-*oo  in  equations  30  and  31.  In  such  case  the  wire  temperature  and  the 
interface  temperature  will  be  equal. 


<9,  (t)  =  0,(  t 


_  16  k2Q'  f"  (l-  e~*2|>2t) 

1  a2b7r2h^  x3  hjn  [P2  +  Q2 


where  ^im  =  xJ0(bx) - j-jJj(bx) 


and  Jim  Qt  =  xY0(bx) - rY.(bx) 

H—  i  o  ah2  1 


By  proper  change  of  variables,  equation  32  reduces  to  Ta  given  by  Friedman  6/. 

Equation  29  may  be  integrated  numerically  to  provide  a  general  solution  for  the  interface 
temperature.  Equation  29,  together  with  equation  1,  may  be  used  to  predict  a  more 
accurate  ignition  delay  time. 
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EXAMPLES 


1.  A  0.1  mil  diameter  Nichrome  wire  in  lead  styphnate.  Figure  1  shows  a  comparison  of 
predicted  explosion  delays  with  those  obtained  from  computer  experiments.  The  solution 
from  equation  31  for  constant  current  heating  agrees  very  well  with  the  exact  and 
quadratic  solutions  6/.  Figure  2  shows  the  effect  of  the  contact  resistance  on  the 
predicted  explosion  delays  of  the  same  wire-explosive  system.  For  values  of  H  ,  the 
outer  conductivity  of  the  surface  between  0.1  and  <»  ,  the  change  in  the  explosion  delay  is 
small,  but  decreasing  H  further  results  in  considerable  increase  in  the  explosion  delay. 

2.  Heating  function  i(t)  =  0.323  sin  (3l4560t).  This  corresponds  to  a  0.1  mil  diameter 
Nichrome  wire  in  lead  styphnate,  the  temperature-time  history  calculated  from 
integrating  equation  30  numerically  is  compared  with  the  exact  solution.  The  explosion 
time  from  equations  29  and  1  was  found  to  be  2.75  nsec.  The  explosion  time  found  by 
numerically  integrating  the  full  set  of  differential  equations,  including  chemical  reaction, 
is  2.74  /tsec.  The  solution  presented  in  figure  3  is  in  agreement  with  those  previously 
published. 
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Compartton  of  predicted  explosion  daisy*  with  thoM  obtained  from  com pu tar  experiments,  0.1  mil  diameter 
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_  i  Friedmans' soproxlmste  solution.  " 

equations  (30)  and  (31)  with  H  ■  1E9 


Figure  1.  Comparison  of  predicted  explosion  delays  with  those 
obtained  from  computer  experiments 


Comparison  of  tamparatura-tlma  history  of  surface  of  load  styphnata  dua  to  sinusoidal  currant  In  0.1  mil  dlamatar 
Nlcnroma  wirai_M.axact  solution!  X  «*C.  aquations  (21)  and  <2t)  with  H  ■  1E9. 


Figure  3.  Temperature-time  history  at  interface  of 
explosive  and  wire  due  to  sinusoidal  current 
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2 

h 

2 

H 

i(t) 

J0<> 

JX( ) 
K00 
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M1 

P 

Pi 

q 

Q(t) 

Q*(t) 
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Oi 

r 
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tra<T) 
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Radius  of  wire 
Specific  heat  of  wire 
Function  of  real  time  t 
Function  of  temperature  T 
Function  of  transform  variable  P 
Diffusivity  of  the  explosive  material 
The  outer  conductivity  at  the  interface  (r=b) 
Heating  current  function 
Bessel  function  of  first  kind,  zero  order 
Bessel  function  of  first  kind,  first  order 
Modified  Bessel  function  of  second  kind,  zero  order 
Modified  Bessel  function  of  second  kind,  first  order 
Mass  of  wire/unit  length 
Transform  variable 
Equation  24 
VrP/h2 

Heat  supplied/unit  wire  length/unit  time 

Q(t)/2  ir  b 

Function  of  P 

Equation  25 

Radial  coordinate 

Outside  radius  of  circle  of  integration 
Time 

Time  at  which  the  interface  temperature  equal  T 
Adiabatic  explosion  time 
Interface  temperature 
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Dummy  variable 

vS/h  2 

Bessel  function  of  second  kind,  zero  order 
Bessel  function  of  second  kind,  first  order 

iVI  i  Cj  /  *■  b 
Function  of  r 
Function  of  p 

Inside  radius  of  integration  circle 
Dummy  variable 

Temperature  of  wire  function  of  time 
Temperature  of  explosive  function  of  time  and  radius 
Laplace  transform  of  temperature 
Dummy  variable 


